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Abstract 

Jet noise is directly predicted using large-scale equa- 
tions. The computational domain is extended in order to 
directly capture the radiated field. As in conventional 
large-eddy-simulations, the effect erf the unresolved 
scales on die resolved ones is accounted for. Special at- 
tention is given to boundary treatment to avoid spurious 
modes that can render the computed fluctuations totally 
unacceptable. Results are presented for a supersonic jet 
at Mach number 2.1. 

1. Introduction 

The recent growing interest in turbulence noise is 
largely due to efforts to develop a high-speed civil trans- 
port plane. The success of this new technology is contin- 
gent upon reducing its jet exhaust noise.The jet noise is 
generated by the tune dependent flow fluctuations in the 

near field which are associated with pressure fluctuations 

that pro pagate to the far field producing the radia t ed 
sound. In theory, direct numerical simulation (DNS) 
based on the compressible Navier-Stokes equations pro- 
vide both the flow fluctuations and the acoustic field. 
However, the resolution requirement for high-Reynolds 
nnmhw turbulent flows makes direct numerical simula- 
tion impractical due to current computer limitations. In 
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large-eddy simulations, the full navier-Stokes equations 
are solved, while realizing that the smaller scales of tur- 
bulence can not be resolved for high-Reynolds number 
turbulent flows. Modelling is used to account for the ef- 
fect of these unresolved scales on the resolved ones. It is 
believed that the large scale structure is more efficient 
than the small scale structure in radiating sound [1-5]. 
This indicates that it is appropriate to perform large-eddy 
simulations to accurately capture the large scales of mo- 
tion while modelling the sub-grid scale turbulence. 

The use of large-eddy simulations (LES) as a tool 
for prediction of the jet noise source has been proposed 
by Manlrhadi et aL [6]. The fluctuating sound source in 
the near field was then used to obtain the far field sound 
through the application of LighthiU’s theory. It was 
shown [6] that because of the non-compactness of the 
source and the need to account for the retarded time, ac- 
curate application of Lighthill’s theory requires, in turn, 
prohibitive com puter storage. The present work is con- 
cerned with avoiding such techniques via direct exten- 
sion of the computational domain to capture the radiated 
field. Large-scale equations are solved for both the sound 
and the acoustic fields. The Smagcrinsky ’s model is used 
for simplicity. But further attention is needed for model- 
ling the effect of the smaller scales on the larger ones, 
particularly for the acoustic disturbance. Careful atten- 
tion is given in this paper to the boundary treatment, 
which is the main issue in predicting an accurate acoustic 
field. 

2. floverning Equations 

The flow field of a supersonic jet is governed by the 
compressible Navier-Stokes equations, and can be de- 
composed into filtered and residual fields, namely 
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f = f+r 


a) 


where an overbar denotes the resolved (filtered) field and 
a C‘) denotes the unresolved (subgrid) one. The mean of 
the filtered field is the mean of the total field. Upon sub- 
stituting this splitting in the full Navier-Slokes equa- 
tions, the Altered compressible Navier-Stokes equations 
in cylindrical coordinates takes the form 
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Here Q is the unknown vector, F and G are the fluxes in 
the x and r directions, respectively; S is the source term 
that arises in cylindrical polar coordinates; and k is ther- 
mal conductivity. The total enthalpy is I, the total energy 
is E, and Cy are the viscous stresses. This system of equa- 
tions is coupled with the equation of state for a perfect 
gas. Here, a tilde denotes Favre fitering, 



The unresolved stresses ty appearing in equations (4-6) 
need to be modeled. 


The Smagorinsky’s model [7] was chosen to repre- 
sent the effect of the subgrid-scale turbulence stresses. 

^ij = fc g (8jj/3) — 2pV^^Sjj — (8) 

where kj is the kinetic energy of the residual turbulence. 
The strain rale of the resolved scale is given by 
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The summation S _ is zero for incompressible flow, v R 
is the effective viscosity of the residual field, 

V R = <C,A f ) ' ( 10 ) 

and Af is the filter width given by 

A f =(A x A r ) 1/2 (11) 

For the beat equation, Edison [8] proposed the eddy vis- 
cosity model 

v r3T 

15 = "pf.s; <12 > 

where Pr,is the subgrid-scale turbulent Prandtl number, 
which can be taken as 0-5. Smagorinsky’s constant C, in 
equation (10) is given by Reynolds [9] as 0.23. In the 
presence of mean shear, however, this value was found 
to cause excessive damping of large-scale fluctuations, 
and in his simulation of turbulent channel flow, Deardor- 
ff [10] used C, = 0.1. Piomelli et al. [1 1] found the opti- 
mum value of C, to be 0.1. This value was used herein. 


Although Smagarinsky’s model is one of the sim- 
plest models, it assumes a balance between die subgrid- 
scale energy production and its dissipation. This balance 
may not be true in several situations. Further work is 
needed for developing more appropriate subgrid-scale 
models for jet noise predictions. 

4. Numerical Scheme 


The importance of the dispersion and dissipation of 
a given scheme, used in connection with the computa- 
tional acroacoustics, was highlighted by Hardin [12]. 
Both effects are crucial in computational aeroacousdcs, 
and can render the computed unsteady part of the solu- 
tion completely unacceptable. As such, high-order accu- 
rate schemes are required for problems in computational 
acroacoustics. 


3. Subgrid-Scale Modelling 
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A fourth-order accurate in space, second-order accu- 
rate in time scheme is used, which is an extension of the 
MacConnack scheme by Gottlieb and Turkel [13]. Man- 



jet computations. 


kbadi ct al. [6] used this scheme to study the structure of 
axisymmetric supersonic jet flow. Ragab and Sheen [14], 
and Faro uk, Oran and Kailasanath [15] have also suc- 
cessfully applied this scheme for the study of nonlinear 
instability problems in plane shear layers. Sankar, Reddy 
and Hariharan [16] performed a comparative study of 
various numerical schemes for aeroacoustics applica- 
tions, and found that this scheme offers high spatial ac- 
curacy. In this scheme, the operator is split into two one- 
dimensional operators and applied in a symmetric way to 
avoid biasing of the solution: 

Q n + 2 = L 2x L 2r L lr L lx Q " (13) 

where L represents the one-dimensional operator. Each 
operator consists of a predictor and a corrector steps, and 
each step uses one-sided differencing: 

Predictor. 
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and likewise for the radial direction. The scheme be- 
comes fourth-order accurate in the spatial derivatives 
when alternated with symmetrical variants. Let L] be the 
one-dimensional operator with forward difference in the 
predictor and backward difference in the corrector, then 
L 2 will be the one-dimensional operator with backward 
Hiffrn-nra in the predictor and forward difference in the 
corrector. At the computational boundaries, flux quanti- 
ties outside the boundaries are needed to compute the 

spatial derivatives, and these can be obtained using third- 

order extrapolation based on data from the interior of the 
domain. 


5 Boundary Conditions 

Special attention is given herein to boundary treat- 
ment in order to avoid non-physical oscill a tio ns dial can 
render the computed acoustic field unacceptable. Several 
boundary treatments were considered [17,18]. The 
boundary treatment discussed below were found to be 
stable, non-reflecting, and most suitable for the present 


5.1 Inflow Boundary 


At the inflow boundary (x=0), the radial boundary is 
split into hydrodynamic disturbance and radiation re- 
gimes, which are treated differently according to the 
physics involved. 

3.1.1 Hvdrod vnamic Disturbance Regime 


At the inflow boundary, a small disturbance is intro- 
duced. This disturbance is assumed to be mainly hydro- 
dynamic in nature, and is specified from the centerline to 
r/D=2. To a first approximation, the inflow disturbances 
is assumed to be small such that the linear stability theory 
applies. A normal mode decomposition for the distur- 
bance is assumed in the form: 

[u’ V p’ p] = 91 {<& (r) e 1 (QX ~ m) } (16) 

The governing equations in this case can be reduced to 
Oir-Sommerfeld equation, which are solved to obtain the 
complex wave number a as the eigenvalue correspond- 
ing to the frequency to and the radial functions <X>(r) as 
the corresponding eigenfunctions. 

*(r) = [u(r) 9(r) £(r) P(r)] (17) 

To obtain the disturbance solution a mean flow must 
be specified. In the present work, the analytical functions 
proposed by Tam and Burton [19] to fit the experimental 
data of Trout! and McLaughlin [20] were used. The mean 
axial velocity is given by 
U = 1 for r<h 

U = exp [-In (2) for r>h (18) 

where b(x) is the half-width of the annular mixing layer, 
and is fitted to the experimental data. The radius of the 
potential core, h(x), is related to b(x) through the conser- 
vation of momentum: 



Velocities are normalized by the stream wise velocity at 
the jet exit centerline; density by the exit centerline den- 
sity; and distances by the nozzle radius. The static pres- 
sure is uniform in the radial direction. The total 
temperature is assumed uniform, which provides the dis- 
tribution for the static temperature 


T = T + -- — (20) 

o 2pc p 

where Cp is the specific heat under constant pressure. The 
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density is obtained from the equation of state. 

For the supersonic regime, all characteristics travel 
in the flow direction. Thus the primitive variables are 
given at x=0 as outlined above. An extrapolation is need- 
ed for the two points outside the domain (x<0). Rather 
than using the conventional polynomial extrapolation, 
we choose the extrapolation function to fit the linear 
growth rate, namely: 

Bx i(a,x-cat), 

u (x, r) «= Ae 9t{u(r)e } (21) 

where P represents the growth rate that can be obtained 
either from the linear theory or from the interior solution. 

In the subsonic regime, the following three charac- 
teristics are specified according to the linear stability so- 
lution: 

p t +pcu t = Cj 


Webb [21]. The pressure condition is the same as that ob- 
tained by Bayiiss and Turkel [22], Enquisl and Majda 
[23], and Hariharan and Hagstrom [24], namely: 

P t = ~ r(9) [r P x + R P r + r] ^ 

However, for updating the rest of the primitive variables, 
Tam and Webb have shown that the momentum and con- 
tinuity equations should be used to account for the pres- 
ence of entropy and vorticity waves at the outflow 
boundary. The spatial differencing used in the inner code 
is employed to evaluate the derivatives which appear in 
equation (26). 


For the outflow regime of large radius with the local 
Mach number less than 0.01, the outflow condition is re- 
placed by the conventional radiation condition of Section 
5.12. 


S3 Outer Radial Boundary 


P t -C 2 p - C 2 (22) 

pcv t = C 3 

The fourth characteristic is outgoing and is obtained 
from the interior solution: 

p t -pcu t ■ C 4 (23) 

The four characteristic equations are then solved togeth- 
er to obtain the time derivatives of the variables, which 
are used to update the solution at the inflow boundary. 

5.1.2 Radiation regime 

In die radiation regime (r/D>2), the conventional 
acoustic radiation condition applies: 

Q t = -r( 6 )[|Q, + I Q r + 2] 

where: 

Q - [o v p p] 

r ( e) 

and M is the local Mach number, c is the sonic velocity. 
The spatial derivatives which appear in equation (24) are 
evaluated in an identical manner as the inner flow deriv- 
atives. 

5-2 Outflow Boundary 

The outflow treatment is based on the asymptotic 
analysis of the linearized equations as given by Tam and 


At the outer radial boundary (r=r — , Ootoq^, 
the conventional radiation boundary condition of Section 
5.1.2. is used. 

5.4 Centerline Boundary 

The results presented herein are for axisymmetric 
disturba nce , for which the boundary condition at r=^) can 
be stated as 

^[u p p] = 0 (27) 

V = 0 

The cc atexiinc treatment for non-axisymmetric distur- 
bances is not obvious, and is addressed in a separate pa- 
per by Shih et al [25] 

6. Results and Discussion 

The numerical results of the flow and acoustic fields 
of a cold, nearly perfectly expanded axisymmetric super- 
sonic jet of Mach number 2. 1 will be presented. The total 

temperature of the jet is 294° K, and the jet exit pressure 
is 0.0515 atm. The Reynolds number based on exit con- 
ditions is approximately 70000. In the present calcula- 
tion. the jet is excited at a Strouhal number of 0.2 with 
the Strouhal number defined as St=fD/U 0 , where D is the 
nozzle exit diameter and U 0 is the jet exit centerline ve- 
locity. The computed results will be compared with the 
analytical solutions of Tam and Burton [19] and the ex- 
perimental measurements of Troutt and McLaughlin 
[ 20 ]. 


(24) 


(25) 
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The computational domain for this problem extends 
axially from x/D=2.5 to x/D=35, and radially from cen- 
terline to r/D=32, as shown in figure 1. Due to the steep 
mean flow gradient encountered at the jet exit, the com- 
putational grid was begun at an axial distance x/D=2.5 
from the actual jet exit The appropriate boundary treat- 
ments are also indicated on this figure. Hie computation- 
al grid consists of 391 equally spaced points in the axial 
direction. In the radial direction, ISO points are used and 
stretched between centerline and r/D=2.5 with concen- 
tration of grid points around r/D=0.5. Between r/D=2.5 
and r/D=16, 130 equally spaced points are used with a 
s pacin g equal to that of the last stretched points. 

Figure 2 shows the instantaneous distribution of p, 
p, u and v. As one can see that the solution is clean from 
boundary reflections. The wave-like nature of the flow 
field is evident The sound pattern seems to peak at angle 
jnf Wnr<»/t by the stream wise position where the fluctua- 
tion reaches a maximum. The spreading of the jet can be 
seen in the instantaneous axial velocity and density dis- 
tributions. 

Figure 3 shows the mean flow Mach number pro- 
files at various stream wise positions. The computed 
Mach number profiles i n d icat e the spreading of the jet, 
but it is underpredicted when compared with the mea- 
surements. This might be attributed to the fact that the 
present c al c ulatio n is axisymmetric, and the uncertainty 
of the subgrid-scale model used as was discussed in sec- 
tion 3. Secondly, the present inflow boundary treatment 
might need to be modified to allow for specifying en- 
trainment. 

Figure 4 shows contours of the root-mean-square 
values of the axial momentum and pressure distributions. 
The preferred emission is evident in the figure. 

Figures 5 and 6 show the axial development of the 
aTiai mom entum and pressure fluctuations on the nozzle 
Up line. The computed results are compared with the 
measurements of Troutt and McL au gh lin [20]. As can be 
seen in figure, the computed curve peaks at about 5 

diameters downstream of the measured peak. It implies 
that the noise source of the calcula t ed results is located 
somewhat downstream of the measured peak. In the ex- 
periment of Troutt and McLaughlin [20], instability 
waves were excited by a single point glow discharge 
pi/wnt^/t flush near the nozzle exit. However, as pointed 
out by Troutt and McLaughlin, the measurements 
showed that the excited motion of the jet is dominated by 
the first helical mode in addition to the axisymmetric 
mode It is expected that the peak of the curve would 
move upstream if the axisymmetric and the helical 


modes excitations are applied simultaneously. 

Figures 7 and 8 shows radial profiles of the axial 
momentum and pressure fluctuations in the jet for sever- 
al downstream positions. The radial peak of the fluctua- 
tion a mplitude moves toward the jet centerline, and the 
width of the shear layer increases with downstream loca- 
tions. The profiles at x/D=5 in both figures are typical of 
the eigenfunctions of linear stability analysis. 

Figures 9 and 10 show the frequency spectra of axial 
velocity and pressure at several downstream locations in 
the ceata of shear layer. The spectra indicates that the 
dominant frequency is that of the fundamental forcing 
frequency and the subsequent harmonics. The ampli- 
tudes of the fundamental and harmonics decay at down- 
stream location, but not so rapidly as observed in the 
measurements. 

Figure 11 shows the sound pressure level distribu- 
tion in the far field for the present calculation, the axi- 
symmetric mode of Tam , and Burton [19], and the 
experimental measurements of Troutt and McLaughlin 
[20]. It is seen that the patterns of the sound pressure lev- 
el c on to urs are all qualitatively similar. The lobed ap- 
pearance of figure 11a indicates that the present 

calculation predicts a strong noise radiation in a direction 
at approximately 30° from the location of the maximum 
fluctuation. However, the present results show a down- 
stream shift of the lobes when compared to the analytical 
entnrinns and experimental measurements. As pointed 
out in the discussion of figure S, the c a lculated distur- 
bance peaks downstream of the measured data and the 
spreading is underpredicted. This implies that the field 
«hap<» of the sound pressure level would be 

di sp lac ed downstream relative to Troutt and McLaugh- 
lin's observations. 

Figure 12 shows the calculated and measured sound 
field directivity at a circle of radius 24D with center at 
the jet exit centerline. The angle is measured from the jet 

exit centeriine. The calculated peak occurs around IS de- 
grees, which is smaller than the measured one. This is 
due to die downstream shifting of the field pattern of the 
overall sound pressure level, as was discussed in figure 
11 . 

7. Conclusions 

Direct prediction of jet noise seems to be feasible for 
the axisymmetric case. In the large-scale simulation pre- 
sented herein, the computational domain is extended to 
capture the radiated field, while modelling the effect of 
the onaiw scales on the larger ones. However, careful 
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attention is needed as was done herein, to boundary treat- 
ment, which becomes a serious issue in computing the 
acoustic field. The only limitation of large-scale simula- 
tion as compared to direct numerical simulation is that 
the smaller scale is not predicted. However, it is believed 
that the larger scales are more effective than the smaller 
sc ales in e mitting sound. In conventional large-eddy 
simulation, the objective is the prediction of the mean 
flow field. But, herein, the objective is the fluctuation 
f ield Thus, careful attention is given to the algorithm 
used and boundary treatment for handling wave-type so- 
lutions. 

Results presented for the axisymmetric field of 
Mad) number 2.1 jet with single frequency disturbance 

at the inflow exhibit features consistent with the observa- 
tions. Three dimensional effects are needed to produce 
jet spreading rate consistent with the observation. The 
small scale modelling might also need crucial attention, 
not only fa adequate spreading of the jet, but also for ex- 
tension to acoustic disturbances. The predicted acoustic 
field, however, seems to be in qualitative agreement with 
the observation. As previously pointed out [19,20], the 
arimiuhal modes seem to be responsible for a higher 
peak emission angle as compared to the axisymmetric 
case. 
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Acoustic Radiation 


Fig. 1 Computational Domain 
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Fig. 2 Instantaneous distributions of (a) density (b) axial velocity (c) pressure (d) 
radial velocity. 



Local Mach ooDbs 


Fig. 3 Mean Mach number profiles at various streamwise positions. 





F»g. 7 Radial profiles of root-mean-square values 


of axial momentum. 


Fig. 8 Radial profiles of root-mean -square val 
ues of pressure. 







Fig. 11a Sound pressure level contours for 
present calculation. 


Fig. lib Sound pressure level contours for 
Troutt and McLaughlin’s experiment 


Fig. lie Sound pressure level contours for Tam 
and Burton’s analytical solution. 
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Fig.12 Directivity of jet noise at R/D=24. 
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